arXiv: 1501.0021 lv3 [hep-th] 16 Sep 2015 


Vanishing beta function curves from the functional renormalisation group 

P. Mati 

MTA-DE Particle Physics Research Group, P.O.Box 51, H-4001 Debrecen, Hungary 
Budapest University of Technology and Economics, H-1111 Budapest, Hungary and 
Eotvos University, H-1117 Budapest, Hungar^ 

In this paper we will discuss the derivation of the so-called vanishing beta function curves which 
can be used to explore the fixed point structure of the theory under consideration. This can be 
applied to the 0{N) symmetric theories, essentially, for arbitrary dimensions (D) and field com¬ 
ponent (N). We will show the restoration of the Mermin-Wagner theorem for theories defined in 
D < 2 and the presence of the Wilson-Fisher fixed point in 2 < D < 4. Triviality is found in 
D > 4. Interestingly, one needs to make an excursion to the complex plane to see the triviality of 
the four-dimensional 0{N) theories. The large-A^ analysis shows a new fixed point candidate in 
4 < D < 6 dimensions which turns out to define an unbounded fixed point potential supporting the 
recent results by R. Percacci and G. P. Vacca in: ’’Are there scaling solutions in the 0(A'^)-models 
for large-AT in L> > 4?” [Phys. Rev. D 90, 107702 (2014)]. 
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I. INTRODUCTION 

There have been many investigations on the phase structure of the scalar 0{N) models p1l2l[8HIO] using functional 
renormalisation group (FRG) techniques [TTl [T^ [TSl HH] in various dimensions. Critical exponents were computed 
with high precision, fixed points were found, and many physical aspects seem to be well understood. Hence, we can 
safely state: the FRG method is a reliable tool for non-perturbative calculations which can be used to consider the 
phase structure of a model and consequently to study e.g. the appearance of spontaneous symmetry breaking (SSB). 
For instance, according to the Mermin-Wagner (MW) theorem a continuous symmetry cannot be broken 

spontaneously below a lower critical dimension, which is D = 2 for the 0{N) symmetric theory for N >2. This has 
been recently shown in the framework of FRG |2H]. Numerical evidence for the MW theorem is presented in |29) . 
The only exception is in H = 2 for V = 2 (the XY model), when an infinite order phase transition is present [3l [4]. 
That is called the Berezinskii-Kosterlitz-Thouless (BTK) phase transition. To detect this kind of symmetry breaking 
the local potential approximation (LPA) is not sensitive enough [^. Throughout this paper we will not consider 
the special case when N = 2 in D = 2 for these aforementioned reasons. Beyond LPA the BTK phase transition 
can be observed in the framework of FRG m\- In summary, we can say that the 0{N) scalar model represents 
an excellent playground for testing new techniques. The goal of this paper is to introduce the so-called vanishing 
beta function (VBF) method in the framework of the 0{N) symmetric scalar field theory, where the Wetterich FRG 
equation mm is used in the LPA assuming that the potential is analytic around the vanishing field. There is an 
extensive study about the radius of convergence of the Taylor-expanded effective potential in the large-V limit of 
the 0{N) model in [22) . Expanding around vanishing field enables us to provide a one-parameter representation of 
the fixed point solution which actually defines the VBF curves through the couplings. This single parameter is the 
quadratic coupling of the theory (m^) and the roots of the VBF curve will define the position of the fixed points at 
a given truncation level. In this approach we will lay down general requirements that a root must satisfy in order to 
define a stable fixed point potential. Using this polynomial expression for the couplings will reveal an interesting root 
structure, which depends highly on dimensionality but there is a dependence on N, too, which is rather quantitative. 
As a consequence of the approximations used in the VBF technique could lead us to potentially artificial hxed points 
even in the higher order of truncation. However, it can be shown that these ’’spurious” fixed point solutions are 
excluded by taking the infinite limit in the order of the expansion, although strictly mathematically they must be 
considered as a valid fixed point at finite truncation. We will show, from a different aspect that has been used in 
[28] and [29], that the Mermin-Wagner theorem is not violated, although at finite truncation order one would draw 
the opposite conclusion. An equally interesting finding is the direct observation of triviality above four Euclidean 
dimensions; however for D = A triviality can be found only if we track down the root structure on the complex plane, 
too. In theories 4 < < 6 a new fixed point candidate is found in the large-V limit, a detailed analysis is shown for 

the D = 5 case. However, it turns out that this hxed point dehnes an unbounded potential from below in agreement 
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with the findings in |35| . 

In the last two sections we will analyse the N dependence and some fractal-dimensional results are presented, too. 
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A. The functional renormalisation group 

The effective average action of an 0{N) symmetric scalar field theory in D Euclidean dimensions using the LPA 
approximation will have the following form: 




d^x 




(1) 


Uk is the dimensionful potential that depends only on the 0(N) invariant terms of the theory, namely on (jr, where 
(j) is the vacuum expectation value (VEV) of the field. The subscript k stands for the RG scale (i.e. the Wilsonian 
cutoff), on which the effective theory is being defined. The dependence of Eq. Q on the RG scale parameter k is 
governed by an exact functional differential equation [mini nulls] 

dtVk = ^Tr (rf ^ + Rk) dtRk. (2) 


Here, we have introduced the logarithmic flow parameter t = ln(A:/A), and a momentum dependent regulating function 
Rk{q^), which ensures that only the fluctuations above the Wilsonian cutoff scale are being integrated out. The A 
is the momentum scale on which our theory was initially defined. T). '[(j)] is a shorthand notation for the second 
functional derivative S'^Tk/S(j)Sip, and the trace denotes an integration over the momentum. Alternatively one can 
introduce a p = and from now on we are going to use this as our variable. Let us note that in the LPA, the 
two-dimensional Wegner-Houghton RG equation (i.e. the sharp cutoff regulator) is mathematically equivalent (see 
e.g. [20]) to the effective average action RG equation nnnaiii] with the power-law regulator Rk{q^) = q‘^{q^/k'^) ^ 
EE] with 6=1 and the functional Gallan-Symanzik RG equation [21j . By inserting Eq. 0 into Eq. 0 one will 
obtain the flow equation for the effective potential 


dtUk = 


1 

2 


d^q 


dtRk 


TV- 1 


q^ + Rk + U'f. 


+ Rfe + C/' + 2pl7" 


( 3 ) 


This equation on the right-hand side defines a loop-integral structure with the propagators of the TV — 1 Goldstone 
modes and the single massive radial mode. The integration by the momentum q can be performed by choosing the 
Rk{q^) regulator function in such a way that it satisfies all the following requirements: P/c approaches the bare action 
in the limit fc —>■ A and the full quantum effective action when fc —>■ 0 |18j . Various types of regulator functions can 
be chosen, but a more general choice is the so-called Compactly Supported Smooth regulator (CSS) [5TH55] which 
recovers all major types of regulators in its appropriate limits. By using a particular normalisation |32L 133] and the 
notation y = q‘^fk‘^, the dimensionless CSS regulator {rk{q^) = Rk{q^)/q^) has the following form 


norm 

CSS 


iy) 


exp[ln(2)c] — 1 


exp 


\n{2)cy^ 
1 — hy^ 


- 1 


0(1 - hy'>), 


( 4 ) 


with the Heaviside step function 9(y) where the limits are 


lim 1 

c— 

lim 1 

c—>-0 


lim 1 

c—^0 


1 

1 


exp[ln(2)y^] — 1 


(5a) 

(5b) 

(5c) 


Thus, the CSS regulator has indeed the property to recover all major types of regulators: the Litim [T5], the power-law 
m and the exponential mKH] ones. Let us note that in the LPA the momentum integral of the RG equation can 
be performed analytically in some cases, e.g. by using the Litim, the sharp, and the 6=1,2 power-law cutoffs. ^ 
using the optimised (or Litim) regulator [11], Rk{q^) = (fc^ — q^) ^(^^ — 9^), one can evaluate the integral in Eq. 0 
analytically. In the present paper we are always going to use this regulator. In the meantime we introduce the 
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following dimensionless variables in order to be able to extract the fixed point structure of the Wetterich equation 
Eq. 11): 


u{p) =U/_k°, 
p = ^ P ^. 


( 6 ) 


At the end of the procedure one will obtain the flow equation for the dimensionless effective potential in D Euclidean 
dimensions: 


dfU = —Du + {D — 2)pu' + {N — 1) 


A 


D 


A 


D 


1 + It' 1 + m' + 2pu" 


An = 


2D-1t^d/2y(^D/2)D' 


(7) 


The hrst two terms arise due to the canonical dimension of the potential and of the fields, whereas the third and 
fourth terms are a consequence of the fluctuations of the N — 1 Goldstone modes and the radial mode, respectively. 
The numerical factor An comes from the angular integration of the T)-dimensional integral in Eq. ([^, and it could 
be absorbed by using a rescaling p —>■ pjAn and u —>■ ujAn, but we will choose not to do that, to show the explicit 
dimension dependence. {An can be set to 1 at any point.) In the following section we are going to study the scaling 
solutions of Eq. Q . 


II. THE VANISHING BETA FUNCTION CURVES 


In the following, we are going to use the most common ansatz for solving Eq. Q. Namely, it is the Taylor expansion 
of the effective potential in the field variable. First we are assuming that we can expand the potential in a Taylor 
series around vanishing held. That is: 


u{p) 



( 8 ) 


For the sake of simplicity we are going to use the following notations for the coefficients \i = uP{0). The scale 
dependence is encoded into these dimensionless couplings. Most of the time we are also going to use for the quadratic 
coupling Ai = rp, sometimes for the quartic coupling A 2 = A and A 3 = r for the sextic coupling. Keeping this in 
mind, we can look at the how equation of the effective potential and differentiate it once, then evaluate it at p = 0 . 
So, we get the how equation for the mass 


dtuP = (D — 2)rP — DrP 


3XAn 
(1 + 


{N-l)XAn 

(f + 


(9) 


If we are looking for the scale independent solutions (i.e. the hxed point solutions) of this partial differential equation, 
one can take dtm? = 0. By doing this, we can express A by using only the mass term 


2m^ (1 + 

{2 + N)An 


( 10 ) 


This curve dehnes the value of A, provided dtm^ = 0; i.e. this relation is only true when the mass stopped running 
US [3D]. If one does not have quartic coupling, i.e. A = 0, then the solutions for this equation are just the roots of 
X{m?)] that is, rp = —1 or rn? = 0. Now we derive the how equation for the quartic coupling, too. To do so, we 
need to perform the same idea as before, but now we need to diherentiate the how equation of the effective potential 
twice with respect to p and only then evaluate it at p = 0. One will have then 


dtX = 2{T)-2)X-DX +An 


ISA^ 


5t 


2 A 2 


(1+m^)^ (l + TO^)^] ^ ^l(l + m^)^ (1 + m^)^ 


( 11 ) 


Again, if we are interested in the hxed point of the equation we need to take dtX = 0, which e nab les us to emress the 
sextic coupling as r = T{X,rP). If we are looking for the hxed points of both equations Eq. ( [ll| ) and Eq. (|^ we can 
express the sextic coupling by only using the mass as an explicit parameter: t = T{X{rP),rrr)= t{tiP). This looks 
like 

2m^ (1 + (D (f + m^) (2 + V) - 4 (2 + iV + 2m^{5 + N))) 

{2 + N)^{A + N)Al 


(12) 
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This function defines the value of r, provided dtm'^ = 0 and dtX = 0. When we are looking for a fixed point for the 
whole system of equations we need to set the Ihs of Eq. (12) to zero (which sets dtT = 0 automatically), thus providing 
the values for w? where the fixed points are found. The general statement is the following: one can express the nth 
coupling simply by using the mass term m? as an explicit parameter. This nested formula reads 


An = A„(An-i(A„_2(...A2(rn^))),rn^) = 


(13) 


It is straightforward to prove this using induction (see the Appendix). As a consequence, one can find a formula 
which tells the general form for the nth coupling for n > 2: 


2 [n/ 2 ] 


n-2 n-2 


An = (-1)"+^ — - m\l + E E 

m-l TT /r,„' , Ar'l[(«-l)/4 


Al-^ n (2* + A^)‘ 
2=1 


z=0 j—0 


(14) 


Here the notation [.] means the integer part of its argument and fij{N) is an integer valued function of where a 
is an integer, too. From Eq. (14) we can conclude that (apart from the prefactor which depends only on N and D) A„ 
is a polynomial of over the integers, which has roots at = 0 and nk = —1 for every n > 2. Another interesting 
consideration is that A 2 is the greatest common divisor of all the A„ VBF polynomials, n > 2. Of course, there are 
many more roots in the complex numbers domain in general (actually the number of roots is growing with n ), but we 
are only going to consider the real ones, for which the physics is meaningful. However, in special cases like in = 4 
and 4 < I? < 6 in the large -N, one can observe a convergent series of the unphysical roots on the complex plane from 


which the physically sensible results can be extracted. It is worthwhile to emphasise that Eq. (14) is the most general 
form of the couplings A„(to^); their qualitative behaviour highly depends on the dimensionality, as we are going to 
see later on. It also depends on the number of the fields N, of course, but this dependence is rather quantitative. 

The coupling A„ = A„(to^) in Eq. (14) defines a curve in terms of w? on which dt\n-i = 0, i.e. the beta function 
of the (n — l)th coupling, vanishes. But we know that A„_i = Xn-iiik) [which we have already in hand, otherwise 
we could not build up A„(m^)] defines a curve on which dtXn -2 = 0, and we can continue this till n = 1. For this 


reason from now on we will refer to the curves defined by Eq. (141 as the VBF curves. 


As a next step we would like to extract the fixed points of the theory. The VBF curves of course define the possible 
values of the couplings for the fixed points. In order to find a fixed point we need to perform the following procedure. 
The VBF A„ defines a curve where dtXn-i vanishes, but it does not say anything about dtXn itself. The curve which 
on dtXn vanishes is defined by A„+i = A„+i(m^), and so on. But we need to cut our Taylor expansion at some order, 
to be able to carry out real computations. Let us say we truncate it at the nth order, but then again we would need 
the beta function of A„ to be zero, which is encoded in the VBF defined by A„+i. Since we expanded the effective 
potential in Taylor series until the nth order, clearly we cannot construct that curve. The only way to get 9tA„ = 0 
is to set A„(rn^) = 0 which gives zero for its beta function automatically. In other words, assume that we would like 
to have an (n — I) order expansion, but before we do so we do not set A„ = 0 for the moment. If A„ were not zero we 
would need its beta function to vanish, too, in order to find a fixed point. But we do not want to compute its beta 
function, because at the end we are satisfied with an n — 1 order expansion. Since we did not set A„ to zero yet, we 
can express it through the vanishing beta function of the (n — l)th coupling because its fixed point equation has the 
form dtXn-i = 0 = F{m^) + A„, where F{m‘^) is a polynomial in m^. Now, A„ = —F{m?), which just defines its VBF 
curve. Since we did not want to have this term (i.e. the nth), we can set F{rn?) = 0, which is satisfied at its roots 
tuqS; hence A„ = 0, too. Thus, although practically we expanded the effective potential till the order n and we were 
able to construct the VBF curves all the way till order n, we must find the roots of the nth VBF curve to set the nth 
coupling to zero, which would imply by construction that 9tA„ = 0. To make a long story short, we have to solve the 
following equation. 


Xniml) = 0 , 


(15) 


where mg represents the roots of the VBF curve A„(m^). We need to evaluate the other n — 1 VBF curves at these 
vn? = mg points to obtain the values of all the dimensionless couplings at the fixed point, and thus define the truncated 
effective fixed point potential. So far so good, but let us suppose we find a fixed point potential which is unbounded 
from below, i.e. one that defines an unstable theory. That obviously must be wrong; hence we need to bring into play 
another restriction: since the potential is a polynomial, the asymptotics (i.e. the boundedness from below for large 
field values) depends on the highest degree term in the polynomial. As a consequence, we need to exclude the m§ 
roots that give A„_i(mg) < 0 coupling, which is the coefficient of the highest degree non-zero term in the polynomial 
expansion. 

In general we can sum up all these requirements in the following way. Let us define the set M = 
{ml e K|A„(mo) = 0 }. A stable fixed point effective potential can be found by substitnting all mg € M into the n — 1 








5 


VBF curves {A„_i(to^), A„_ 2 (m^),A 2 (to^), Ai = m^} provided A„_i(to§) > 0. [If it happens to be zero, too, one 
needs to apply this rule for the VBF A„_ 2 (m^), and so on.] In other words, the set of TOq’s which defines a true fixed 
point is 

M* = {mg e K|A„(too) = 0 a A„_i(mo) > O} . (16) 

In the following, some example are presented individually foriA<2,2<Zl<4 and D > 4, starting with D <2. 


III. VBF CURVES FOR D<2 


A. Continuous symmetries {N > 2) 


The Mermin-Wagner theorem essentially states that no spontaneous breaking of continuous symmetry is present in 
systems of H < 2. In [2^ numerical evidence was given in the framework of the FRG that the MW theorem indeed 
holds; here a beyond-LPA scheme was used, where the wave function renormalisation was taken into account, too. In 
[25] it was shown, using analytical considerations, that the MW theorem is not violated even at the LPA level. We 
will verify the MW theorem for the case of the expanded potential (where this statement cannot be seen directly) by 
applying the rules that have been settled above. Although the title of this section suggests that we include the case 
of a system with D = 2 and N = 2, we will perform our analysis using D = 2 with N > 3, since in the former case an 
infinite order BTK phase transition is present mm), as it was mentioned above. To detect this kind of symmetry 
breaking the LPA is not sensitive enough (cf. [25]). However, beyond LPA the BTK phase transition can be shown 
in the framework of FRG m\- 

Using Eq. (141 for D = 1 and D = 2 case one can find a simplified expression for the nth coupling 


A„ oc (-1) 


n+1 


m 


n-2 

E 

2=0 


(l + m^)"> g,(iV,iA)(m^)L 


(17) 


Here we only indicated the polynomial structure in rn?, although the prefactors are slightly modified, too. The 
coefficient functions gi{N,D) are defined as follows: 

n—2 

g.(V,U) = ^/,,,(iV“)UU (18) 

3=0 

Interestingly, when setting D < 2 the finding is that all the gi’s are positive for every term, at least till the highest 
order (n = 45) that has been considered in the expansion; this is the observation. This fact suggests that the roots 
(which are not complex) must be either negative or zero. 

Let us find the fixed points according to the rule that has been established in the previous section, using N = 3 
and D = 2. The following VBF curves are found: 

A2 = 

A3 = 

A 4 = 

As = 

Gonsidering only the real roots, we can find that they are all situated in the interval [—1,0]. In Fig. one can see 
the plot of the VBF curves versus . The mere fact that we can find roots is against the Mermin-Wagner theorem, 
because for that to hold true, we should have no roots at all, except at the ending points of the interval [—1,0]. 

We found true fixed point potentials at a finite truncation level according to Fig. [^ which should not be there as 
a consequence of MW theorem. The situation is getting worse when we go to higher order in the truncation: indeed, 
in this case we are going to have more zeros; hence more and more fixed point potentials emerge. How can we resolve 
this contradiction? One way to overcome this situation would be that if all the roots between —1 and 0 turned into 
complex valued roots as we went to higher order. 


— ^m^(l -I- m^)^, 
o 

647r^ 9,. 2\3 


175 

5127r3 


7875 


(l + m")'^(5 + 27m"), 

m^(l -k rri^f (25 -k 670m^ + 1671{my 


409671 ^ 

606375 


m^(l -f rri^f (175 -k 18595771^ -k 161115 (to^) 2 -k 254799 (to^)^) , 


(19) 
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FIG. 1: The VBF curves for the couplings A 2 , A 3 , A 4 and A 5 in the 0(3), D = 2 model . Each curve is defined by the 
vanishing beta functions: dtrn^ = dt \2 = dtX^ = dtX^ = 0 respectively. The roots for X^{m?) are indicated by black 
dots. These roots are going to define the fixed point potentials at the truncation level n = 4. Two of its roots are 
indicated by a green dot and a red triangle to demonstrate a valid and a false fixed point, respectively. At the green 
dot {m? « —0.484) we find that A 4 > 0, and evaluating A 4 at the red triangle (m^ = —0.137) we get A 4 < 0, hence 
dehning an unbounded potential (see Fig.[^. On A 4 and on A 5 a scaling is performed for the sake of the presentation. 



FIG. 2: The fixed point potentials evaluated at the hxed point discussed in Fig. [TJ A stable potential can be obtained 
for « —0.484 and an instable one for = —0.137, which are indicated on the figure by rtst, Um, respectively. 


However, from the order of the expansion we used, we cannot put our hope in this. There is another scenario, 
too. Let us analyse the VBF curves more carefully. Apparently, one can find a pattern of the roots for the VBF 
curves: evaluating the A„_i VBF at each of the roots of A„, one will find real numbers alternating in sign starting 
from the closest root to —1 with sgnA„_i(mi) = +1 (here the 1 in the subscript indicates the closest root of A„ to 
— 1). For instance in our example in Fig. that was « —0.484 corresponding to the stable potential (Fig. [^. 
In other words, from the set M defined above, we will consider only M \ {—1,0} and the claim is as follows: if we 
make an ordering from the smallest value to the largest in this set, then every odd element of M \ { — 1 , 0 } is in the 
set of M* \ { —1, 0}, hence for every such element A„_i > 0 thus defining a stable potential. So, we are still going to 
have fixed points even though the number of them is halved by taking into account only the fixed points which define 
stable potential. It seems that this procedure does not help too much, but, actually, by doing this, we can extract 
some useful information: each root of the A„ VBF is surrounded by the roots of the A„_i; otherwise the alternating 
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signs that were explained above could not be possible. So, it means that by deriving higher and higher order VBF 
curves the position and the number of roots change in the way that all of the previous roots are around the new ones; 
see Fig. Let us call this interesting pattern of the set of roots the M* pattern for future use. We can do one thing 
with this without knowing anything about the structure but the root pattern statistics: we can simulate a sequence 
of sets of points which behave in this way. 



FIG. 3: The roots of the VBF curves A 2 ,..., A 15 . Observe the pattern which is given by the following rule: any root 
of A„ is between the roots of A„_i (except for —1 and 0 of course). 


Let us consider a randomly generated number Xi which can take a value in the interval (—1,0). We generate this 
number and then we consider two new random numbers: X 2 and V|. The first one can take any value in the interval 
(—l,Vi) and the second in (Vi,0). After we generate values for X 2 and A|, they are going to be the new ending 
points of the intervals where we define again random numbers but this time three: Ag, A| and Af. We continue this 
procedure with the random numbers A^ with i = 1 , 2 , ...,n and n —>■ 00, where n indicates that we defined them in 
the nth step in the interval {X^\, X^\), including —1 and 0, too. By increasing n we can obtain the distribution 
of the points created in the way described above; this is shown in Fig. From this construction one can see that 
the distribution of the randomly generated points is changing in a way that they are accumulating at the two ending 
points of the interval. This suggests that the limit of the probability density for A„ (at least in distributional sense) 
is 


lim fx^{x) 

n—>-oo 


i((5(x + l) + 5(a:)). 


( 20 ) 


The aim was to simulate the zeros of the VBF polynomials A„, and we found that the roots are accumulating at —1 
and 0. Thus, if there exists a limit for the VBFs (lim„_>oo A„ = A), then this will give zero only at —1 and 0, even 
if this limit is not a continuous function. But this is physically well justified since this would indicate that we do 
not find any roots other than the two ending points of the interval, which leads us to the conclusion that no SSB is 
present in dimensions D = 2 with A = 3. Indeed, this is a seemingly paradoxical result, namely that we would expect 
infinitely many fixed points as the degree of the VBF polynomials grows, but at the end we will find two actual fixed 
points defining only one phase, in agreement with the Mermin-Wagner theorem that we wanted to show. We can 
make another remark on the role of the convexity (or IR) fixed point = — 1. From the VBFs one can see that 
every higher coupling vanishes at this point, i.e. A2(—1 ) = A3(—1 ) = ... = 0 ; hence the potential, which is defined by 
this fixed point, is an unstable one, corresponding to the dimensionless potential u = —p^. One can now speculate 
whether this fixed point is a real one or if it is just an artifact of the approximations that we use in the FRG. Similar 
results can be obtained for the D = 1 (see Fig. [^. For N dependence and fractional dimensions, see Secs. VI and 
Eg respectively. 
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FIG. 4: Histograms of the distribution of the points generated in X„: (a) the distribution of Xiq, (b) the 
distribution of X 50 , and (c) the distribution of X 28 ooo- Note that the positions of the points are tending to —I and 0 
and at very high n the number of points between the two end points is negligible; in other words, the probability of 
finding a point between the two end points is tending to zero as n —>■ 00. 



FIG. 5: Roots oi a D = I-dimensional theory with N = 3 field components. A similar structure can be observed as 

in Fig. ^ 


B. Z 2 symmetry {N — 1) 

In the case of the discrete symmetry the 0{N) theory is equivalent to the Ising model, and as the symmetry is 
being non-continuous, the MW theorem does not necessarily hold for = 1. Indeed, it is well known that in the 
Ising model we can find a spontaneous symmetry breaking even in the D = 2 case, which was carried out in an 
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exact calculation by Onsager [5^. Here, we are going to see that we can reproduce this result using the technique 
introduced above. 




FIG. 6 : In the upper (lower) panel the position of the roots of the VBF curves are shown for iV = 1 in £) = 1 
{D = 2). In the one-dimensional case one finds the same pattern that was shown above for iV = 3 with a 
convergence to = — I. In two dimensions the roots are converging towards = —0.835 signaling a SSB phase. 


In Fig. I^we can see the position of the roots for D = 1 and = 2 in a separate plot. The VBF polynomials 
are described by Eq. (17) with different gi coefficients, of course. We can see qualitative difference between the two 
figures; namely, in the one-dimensional case, the first root is converging to = — 1 and for the two-dimensional case 
the asymptotic goes to = —0.835 ± 0.001, defining a threshold for all the roots in higher orders. One can fit a 
function of the form of f{x) = a + on the positions of the first roots, with fit parameters a = — 1.00 ± 0 . 0002 , 
b = 2.08 ± 0.002 and c = 0.85 ± 0.0006 in £> = 1; and a = -0.835 ± 0.001, b = 1.56 ± 0.004 and c = 0.738 ± 0.003 
in D = 2. Now, since the root position pattern is the same as we observed above, we can apply the random number 
procedure for H = 1 in the interval [—1,0] and for = 2 in the interval [—0.835,0]. The limiting result is the 
following: in one dimension the roots are accumulating at m? = — 1 and = 0 ; thus we can conclude only the stable 
Gaussian fixed point potential is well defined at the latter. In the case of the two dimensions we will find two fixed 
points, the Gaussian one {rn? = 0) and another one which is at m? = —0.835. Now, it is a question whether this 
fixed point is stable or not. The finding is the following: for every truncation in the present computation (where the 
highest order was n = 45) Xn{m^ = —0.835) > 0, thus we can safely state that it defines a stable fixed point potential, 
hence indicating a spontaneous symmetry breaking phase. It is interesting that even at LPA level one can observe 
the SSB for the two-dimensional Ising model contrary to the finding in |28j , where even with the numerical spike plot 
technique mui] it is undetectable, due to its dependence on the numerical precision. 
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IV. VBF CURVES FOR 2 < U< 4 


The 0{N) models which belong to the class 2 < U < 4 have the richest fixed point structure; hence studying them 
is the most challenging. We are going to perform a detailed analysis for the only integer dimension found in this 
interval, D — 3 restricting ourselves to iV = 2. We expect here to obtain the well-known Wilson-Fisher fixed point 
in the broken regime, however, we are going to see that there is no clear root structure to be observ ed f or the roots 
> 0. In this case we have the following generic form for the VBFs in 2 < < 4, similarly to Eq. ( (TtI ): 

n — 2 

A„ oc (-l)"+i(m")i+®("-4)(l + ^ g,{N, (21) 

i=0 


Here, gi{N,D) is defined in a similar way as above in Eq. (18), but now, as we can see, the exponent of has 
changed, hence in the sum Eq. (18) there must be rather than , where a{j) represents the remaining 

exponents, after the prefactor of the sum is factored out. The 9 function in the exponent is just the Heaviside step 
function: 0(n — 4) = 1 if n > 4 and 0 otherwise. The first few polynomials from the VBEs for D = 3 and N = 2 are 


A 2 = —37r^TO^(m^ -I-1)^, 

A3 = 37r^m^(TO^ + 1)^(11 to^ + 1), 

A 4 = + lf{23m^ + 4), 

As = -I- 1)® (2993 (to^)^ + 719m^ + 14) , 

Ag =- + 1)® (97167(to2)3 -p 27418(771^)2 + 997771^ - 14) , 


( 22 ) 


Contrary to the cases of H < 2, the coefficients gt can take now negative values, too. Erom these considerations one 
can already expect a different root structure from that which we had for D <2. In Eig. [^we see the VBF curves up 
to XQ^rri^) and in Fig. [^the root structure up to order n = 41. For our analysis we are going to separate the real line 
of m? into three regions. First, we will consider the roots in the interval [—1,0]; then we will turn to the complement 
set, separately for the intervals (— 00 ,—1) and (0,oo). 



FIG. 7: In this plot the VBFs up to Ag are shown for U = 3, iV = 2. For each Ai the root that is closest to —1 (from 
above) is indicated with a black dot. Note that if we consider the position of these points as a sequence then they 


will converge to a finite value = —0.23; see also Fig. 
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From Fig. 10 it is clear that in the interval [—1,0] we have a very similar pattern for the roots (above some n) to the 
one that we had in the two-dimensional case (see Fig. |^. However, there is a striking difference between the D <2 
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case and the present Z1 = 3 theory: the roots seem to be accumulating around the value « —0.23, rather than —1. 
Although this is of course an approximated value, for the sake of simplicity, we will use the equal sign in the following 
whenever we consider this value: m? = —0.23. If one restricts the pattern to the [—0.23, 0] interval, then the random 
number generating model for the root pattern statistics becomes available again. Since the distribution of the points 
signals the accumulation at = 0 and = —0.23 the probability density of finding a root in this interval will have 
the same form as in Eq. ( |20l ); the only difference is that now the Dirac deltas are centred at = —0.23 and = 0 
rather than = — 1 and = 0. We can also check the convergence of the roots by establishing an ordering among 
them: ml, ml and m§, where the first indicates the root that is closest, the second is the second closest and the third 
is the third closest root to = — 1 at each order. The root labeled by ml starts to go towards —1 but at order n = 6 
it stops and starts to converge to m? = —0.23 in an oscillatory way. 



FIG. 8: Here we present the root structure of the VBF curves up to order n = 41, in the 0(3), D = 3 model. 
Observe that we can distinguish three regions where the roots are positioned: > 0, m^ < —1 and G [—1, 0]. 

For a detailed picture see Fig. 
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FIG. 9: The absence of the Wilson-Fisher fixed point at the truncation order n = 8. We can see that the roots of Ag 
to the left define a negative valued \%{mlyp), which will provide an unphysical potential. Moreover, the second root 
to the right does not give a sensible potential either. This phenomenon occurs because of the oscillatory nature of 
ml, but at n —>■ oo we expect to see a stable fixed point potential at the WF fixed point. As has been rescaled with a 

factor of 2 X 10“^. 


On the other hand, the roots ml, ml also converge to this value, and it is possible to fit a curve on them, but this 
time the fit is an exponentially decaying one, as it can be seen from Fig. [TOa] The convergence of these roots confirms 
what we expect from the M* pattern: the two roots that will remain in this interval for n —>■ oo are m? = 0 and 
m? = —0.23, which correspond to the Gaussian and the Wilson-Fisher fixed point, respectively. One needs to check 
whether the WF fixed point defines a bounded potential. 
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n n 

(a) The region —1 < < 0. (b) The region m? < —1. 



n 

(c) The region > 0 


FIG. 10: The root structure of the D — 3, N = 2 theory separated into different regions. In (a) and (b) we can 
conclude for the n —>■ oo limit; that is, we will have the following fixed points: Gaussian, Wilson-Fisher and the IR. 
These suggestions are also confirmed by functions. Regarding (c) for > 0, no clear answer is fonnd: the red lane 

indicates the last real root for m? > 0. For details see the text. 


Substitnting the value = —0.23 into the VBF polynomials A„ will give the following resnlt: there are VBF curves 
for which A„(to^ = —0.23) < 0 and for which A„(m^ = —0.23) > 0. This is due to the oscillatory behaviour of the root 
nil, but we know that this will converge to a finite value, and finally the WF fixed point can be defined as the limit 
of nil- But until this happens, we will find situations which would give the weird result that the Wilson-Fisher fixed 
point defines an unbounded potential from below; hence, strictly speaking, we should not consider it as a physical 
fixed point. This is the case for example with Ag and Ag (see Fig. |^; hence this tells us that there is no Wilson-Fisher 
fixed point at the truncation level n = 8. Even though this is true, we must not forget that the Taylor expansion is 
an approximation, so the real fixed point structure of the theory will be found when n —>■ oo; hence the ’’absence” 
of the WF fixed point can be considered as an artifact of the expansion. Let us consider now the region < — 1. 
A magnified picture of it can be seen in Fig. |10b| Here, we can observe roots which are running into the convexity 
fixed point as n grows. The position of the roots would show a similar pattern to the M* pattern; however, here we 
have difficulties with the unbounded interval (— oo, —I]. We will show how it is possible to overcome such a situation 
in the next section; however, there will be some requirements, which are not fulfilled in the present case. Here, we 
can only assume that all such roots will converge to w? = — I providing us the IR (or convexity) fixed point in the 
n ^ oo limit. For two such roots we can show the convergence by doing a fit. Interestingly, in this case the roots 
are following a different trend to approach —1: the first one is an exponential function f{x) = a + with 

a = —1.014 ± 0.005, b = —0.39 ± 0.01 and c = —1.43 ± 0.05; the second is power-law g\x) = a -I- c(x — 20)^ with 
a = —0.885 ± 0.003, b = 1.35 ± 0.08 and c = —7.2 ± 0.06. Neither of them goes to —1 precisely; however, these are 
the best fits that can be found. For large n values, when these roots are very close to —1, they develop an imaginary 
part; hence they cannot be considered as true fixed points, strictly speaking. One has to go so close to —1 that this 
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effect can be considered as negligible, taking into account that the VBF curves around —1 are extremely flat thanks 
to the (1 + factor in Eq. (14). Thus, finding roots around —1 is not always a reliable thing; it might depend on 
the precision of the root finding algorithm, too. 

The remaining region that we need to consider is the half interval > 0. The position of the roots for this region is 
shown in Fig. 10c One can find again some pattern which could remind us of the M* pattern; however, in this case 
besides the fact that the interval is unbounded from above, the roots do not heave a clear bound even from below. 
The general behaviour of the root positions is that they have the last real value at about m? = 0.01 — 0.02 (indicated 
by the red lane in Fig. 10c), and below that they will have complex values. Even considering the generated complex 
roots does not help us to understand this part of the root structure, however, in the case of the four-dimensional 


0{N) models we can use these complex roots to capture the real physics, as we will see in Sec. VB 


The most important result of this section is the appearance of the Wilson-Fisher fixed point. It may well be possible 
to find more physical fixed points in other fractal dimensions between two and four dimensions for various N values, 
but it is beyond the scope of the present study. 


V. VBF CURVES FOR U > 4 


So far we considered D = 1,2 and 3. In this section we are going to investigate the theories in U > 4. Let us 
look at the first two VBF curves (A 2 and A 3 ) in arbitrary dimensions and field components. We already derived their 
formulas in Sec. |nj Eq. (10) and Eq. (12). Now, we are interested in the roots of the general expressions: 


2 m^ (1 + 

(2 -|- n)A]j 

2m^ (1 + (D {l + m^) {2 + N) - 4 {2 + N + 2m^{5 + N))) 
{2 + N)^{4 + N)Al 


(23) 


for A 2 and A 3 , respectively. The first equation gives zero at = — 1 and = 0 independently from the dimension. 
Among the roots of the second equation in Eq. (23) of course we discover again = — 1 and rri^ = 0, but factorising 
with respect to these roots, one will be left with the equation for the third root, and it will depend on D. Let us solve 
it for TO^. The result is the following: 


-DN -2D + 4N + 8 
DN + 2D -8N -AQ' 


(24) 


This expression shows the ’’running” of the root (the Wilson-Fisher fixed point) with the dimension. For substi¬ 
tuting D = 1,2 and 3 (and N = 3), one would get the results that we obtained in the previous sections. However, 
evaluating it at U = 4, we will lose the N dependence completely. In other words: no matter what 0{N) model we 
consider, the WF fixed point coincides with the Gaussian in the D = A case. In Fig. [TT]one can see the A3 VBF for 
various dimensions. The observation is the following: the only root in the interval (—1, 0) drifts towards the positive 
real values. At 7A = 4 it merges into zero; hence, besides the convexity and the Gaussian fixed point, there are no 
further fixed points present. Above four dimensions the root gets positive values. We can address the question if 
this is a true fixed point or not. If one substitutes the value of these positive roots into the VBF one will 

immediately see that A 2 (mQ) < 0, thus defining an unstable potential at this truncation level. As a consequence, there 
are no other true fixed points besides the Gaussian and the IR, hence we can call the system trivial. It is possible to 
derive the dimension dependence for other Ai’s, too, with similar qualitative behaviour: the roots start from inside 
the interval (—1, 0) with D = 1 and as the dimension grows, they are tending out from the interval. As D takes the 
value of 4, the last (i.e. the closest to = ~1) root of the VBF under consideration merges into rn? = 0, and all 
the others have already obtained positive values; see Fig. where we illustrate this situation with A 4 . 

In the following, we are going to discuss the VBF curves at a higher truncation level for U = 5 in detail, as a repre¬ 
sentative of theories for D > A. However, we need to distinguish the case when D = A, since showing the triviality in 
this case is not that straightforward as it is for D > A-, hence, the four-dimensional case is going to be presented in a 
separate section. In the last two sections various examples for theories in dimensions D > A are shown. 
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FIG. 11: In this figure one can see how the A 3 VBF curves change with the dimension. Observe that the only root 
from the (—1,0) interval runs into zero as H —>■ 4. For D > A the root obtains a value mg > 0; hence there is no 
fixed point in the SSB phase but there is one in the symmetric phase, which turns out to be unstable because 
A 2 (™o) ^ 0. Note that we indicated only integer values for D but the transformation is continuous in the dimension; 
hence all the values between the integers would define similar curves. For this illustration N was set to 3 and the 
VBF curves are rescaled as Ag (x AglO”", where a = 0,4,4, 5 and 6 for H = 1,2, 3,4 and 5, respectively. 



FIG. 12: In this figure one can see the different A 4 VBF curves for different dimensions. Observe that the two roots 
from the interval (—1,0) runs to the right as D grows. For D = A the root indicated with a triangle has got already 
a positive value and the one which was indicated by a dot has just melted into zero. For D = 5 both of the roots 
have positive values. The triangle is not present in the plot, since it took the value m? = 0.673. For this illustration 
N was set to 3 and the VBF curves are rescaled as A 4 oc A 4 l 0 ““, where a = 5, 6 , 7, 8 and 9 for U = I, 2, 3,4 and 5, 

respectively. 


A. Triviality of the 0{N) model in D > 4 

For any D > A theories the finding is that the VBF curves have the following structure: 

n—2 

An OC (-l)”+^(m^)(l + m^)” ^ g,{N, D){rri^y. 

2 = 0 


( 25 ) 
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Notice that it has essentially the same structure that we had in the cases D <2, but there is one crucial difference. 
For theories D < 2 we found the coefficient functions g^’s to be always positive, hence defining only negative (and 
complex) roots for the polynomial; moreover they were inside the interval (—1,0). In this case likewise to 2 < II < 4 
we can obtain negative values for g^’s, too, but contrary to that case, here we do not obtain roots inside the interval 

(— 1 , 0 ): all the real roots in this case are positioned in the disjoint union of the complement set of (— 1 , 0 ), i.e. in 

(—oo, —1] U [0,oo). The first few VBFs are the following: 

A 2 = —2ATT^rn?{m^ + 1)^, 

^3 = m^(TO^ + l)^(39m^ — 5), 

A 4 = — to^(to^ + I)"* (576(m^)^ — 425m^ + 25) , 

As = TO^(m^ + 1)^ (10233 (to2)3 - 176625(771^)2 + 36125to2 - 1225) , 

539 

(26) 



FIG. 13: The VBF curves A' (i=2,...5,8). Observe that the roots of each curve are always outside the interval 
(—1,0). For As we have negative roots, too. For display purposes the VBF curves are rescaled as A' cx Ail0““, 
where a = 1, 3, 6,9 and 19 for i = 2,3,4, 5 and 8, respectively. 


These curves are shown in Fig. |13| Now, we are going to analyse their root structure. In the case for theories in 
H > 4 we can clearly identify a pattern of the roots again, just like we did it for D < 2, i.e. the M* pattern. For the 
particular case D = 5, N = 3 one can see the position of the VBF roots in Fig. What we can observe is just the 
reverse of what happened in the 0{N) models for D < 2. The roots are situated only outside the interval (—1,0), 
and they have a similar pattern to the one we had for the D < 2 cases; i.e. each root of A„_|_i is surrounded by the 
roots of A„. 

We can call it only similar, since we clearly have a problem in the present case: we are not able to use our 
random number generating model that we did in Section [m] for simulating the positions of the roots, because the set 
(— 00 , —I] U [0, 00 ) is unbounded. For that reason, the highest (lowest) roots of the VBF A„+i do not have an upper 
(lower) neighbour from the set of the roots of the VBF A„. However, we are able to do the following trick: let us 
one-point compactify the real line M, meaning that we ’’glue” together the points m? —)■ ± 00 , thus creating a closed 
set for [0,-1]. This may look like an ad hoc idea but we will see that it works. Let us consider for instance the 
orders n = 6 and n = 7 in Fig. 14 The highest valued root in the region vn? > 0 from n = 7 would need a root from 


71 = 6 to restore the right pattern, but apparently there is no such root. However, if we extend our picture with the 
compactified real line, we will see that in the region rri^ < — 1 the right root appears just where we would need it at 
order 71 = 6, and thus the pattern goes on. This is true for all the other roots (see Fig. 14), thus the compactification 
of the real line is well justified in this way, hence providing us with an M* pattern on the compactified real line. 
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FIG. 14: In the upper (lower) panel the roots of the VBF curves are shown for A 3 ,Aio (Ae,A 15 ), for > 0 
{m? < —1). One can observe a pattern of the roots similar to the D <2 case: each root of VBF A„+i is surrounded 
by the roots of A„. The only problem in this case is, since the interval is unbounded, this pattern is not entirely 
true: indeed, the highest (lowest) roots of A„+i now do not have an upper (lower) neighbour from the roots of A„. 
We can overcome this difficulty by one-point compactifying the real line K. For details see the text. 





FIG. 15: Histograms of the distribution of the points generated in V„: (a) the distribution of X 20 , (b) the 
distribution of Xioo, and (c) the distribution of Visooo- Observe that the points are accumulating at —1 and 0 but 

this time from the outside region. 


Ill The result of such a random 


In the interval [0,-1] one can now use the technique that we introduced in Sec. 
number procedure is shown in Fig. |15| 

We can see that the positions of the roots are again accumulating at —I and 0, but now they are approaching from 
the complement interval of (—1, 0). In this way we can see how the triviality is emerging in the limit n —>■ 00 . 


B. Triviality of the 0{N) model in D = 4 

In this section we are going to discuss the D = A 0{N) models, and special attention will be given to the case 
N = 1, which has been in the centre of interest since triviality is predicted to occur in the theory. To have a clue 
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on triviality beyond PT we still need to rely on lattice simulations, which actually support the trivial behaviour of 
such theories [34]. Here, we are going to show a result which suggests that indeed the trivial scenario holds for such 
models; i.e. no UV fixed point different from the Gaussian is present in the 0{N) models taking into account all the 
symmetry respecting terms in four dimensions (the (j)'^ case is shown in Fig. jllj ). We are going to present the result 
for iV = 1 but this holds for general N. The VBF polynomials in I? = 4 have the form 

n-2 

A„ oc + m^)" ^ g,{N, D = 4)(to2)*-®(”-3). (27) 

2 = 0 


In particular the first few for = 1 are 


A 2 

A 3 


A 4 


A 5 


647r^ 

81927r4 


m^{m^ + 1)^, 

+ 1 )^, 


5242887r6 
35 

671088647r® 

315 


{m'^) (m^ + 1)^(14 to"^ — 1), 

(to^)^(77i^ + 1)^ (l68(m^)^ — 41m^ + l) , 


(28) 


By analysing the root structure of these equations one will find that all the roots are situated outside the interval of 
(—1,0) likewise in the H > 4 cases; the only difference is that at the order n = 7 the VBF curve Ay develops roots 
on the complex plane, too, with Re m? > 0. We have already met such a situation in theories defined in 2 < H < 4 
dimensions. In that case, we neglected these roots since we considered them to be unphysical. In the present case 
we still stick to our convention; that is, we neglect them as possible fixed point candidates. However, we can identify 
an interesting behaviour for such nonreal roots. Let us look at the root structure of the theory in Fig. |17[ where we 
present the position of the real roots for all relevant intervals of m?. 




FIG. 16: The real roots of the four-dimensional 0{N) model are shown. No root in the interval (—1, 0) is found; on 
the contrary, all roots are outside this interval. One can see that, just like in the case of U > 4, the roots are 
approaching —I and 0 from the outside region. However, more careful analysis is needed in this case. For details see 

the text. 

What we can observe is that the roots are approaching the points —1 and 0, just like in the case for D > 4. 
However, as we indicated above, complex roots are emerging: the first one and its conjugate from n = 7, two and 
their conjugates from n = 11, three and the conjugates from n = 15 and from n = 18 four complex roots plus 
conjugates. Among these roots, being complex, we cannot make an ordering; however we are able to do that for the 
real part of them. The real parts of the roots are indicated in Fig. Gomplex roots with negative real parts do not 
occur; hence it is enough to consider the Re m? > 0. 

The hgure shows that for the real part of the roots, there is indeed a pattern, namely almost the same that we 
observed for both the cases D < 2 and D > A; the only differences are in the orders where the complex roots 
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FIG. 17: The real parts of the roots. The black dashed line indicates the average threshold between the purely real 
and complex root values. It is at about = 0.022. For every value of Re under this line we have an imaginary 
part for the root. New complex roots are emerging at orders n = 7, n = 11, n = 15 and n = 18. At these orders we 
can see a ’’violation of the pattern” which is being indicated by the arrows: pointed black arrows for the purely real 
and dashed gold for the complex roots. Apart from the newly emerging complex roots the pattern holds. For details 

see the text. 


appear. Here, we observe that above and below the real part of the newly appearing complex roots the pattern 
continues; only at those particular points do we find the breaking of the pattern (these are indicated by red circles 
in Fig. [T^: no root from the real part of the A„_|_i goes between the corresponding roots of A„. Now, if we consider 
the limiting distribution of such a pattern of the points, we will find that this anomaly will not have an effect on 
it: after compactification of the real line, as we did it for the D > 4, the points are just accumulating at —1 and 
0; hence we will find exactly the same result as in Eq. (201, but now only for the real part of m^. We are not 
done yet, since we have only considered the real part so far. This result on the limiting position of the roots does 
not make any sense if we find a finite imaginary part of in this limit. Let us consider therefore the roots on the 
complex plane. This is shown in Fig.|18| We can see the developing imaginary part at the order n = 7 for the first time. 




Re Re m7 


FIG. 18: The roots on the complex plane. In the upper panel one can see the emergence and disappearance of the 
first complex root (rn\). It appears at n = 7 and tends to zero as we go to higher orders in n. In the lower panel the 
second and third complex roots are indicated (m^ and ra\). They appear at orders n = 11 and n = 15 respectively, 
and coincide with the curve defined by the complex root m\. The fitted power-law curves are defined in the text, 

and the fitting procedure was performed on 24 data points. 
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Interestingly, the absolute value of the imaginary part will have a maximum at some point and it is tending to zero 
just like the real part. We can also express the imaginary part as a power-law of the real part close to the origin: 
Im = ±a(Re m?)^, where the parameters are a = 0.162 ± 0.001 and b = 0.589 ± 0.001. Now, considering the 
’’running” of the second and the third complex roots which are appearing first at the orders u = 11 and 

n = 18, respectively, we will find that they behave in a very similar way to the first complex root {ml); moreover, 
they collapse onto each other close to the origin (see Fig. [^, thus approaching 0 with the same exponent. We can 
see that both the real and the imaginary parts of the roots are approaching zero as the order n grows, hence we can 
say that the roots, although they are not defined as physical ones when they have complex values, accumulate at —1 
and 0 in the n —> oo limit. It is interesting that only in the limit is the imaginary part absent entirely, and until we 
take this limit we will have a gap between m? = 0 and the last noncomplex valued root in the particular order n. 
With this procedure we were able to show that only two fixed points are present (with the Gaussian as the only stable 
one) in the unexpanded 0{N) model just like in U > 4. Altogether this signals the triviality for the 0{N) theories 
when the dimension is D > 4. Hence, we were able to show a non-perturbative evidence of the triviality for theories 
defined in I? > 4. 

Regarding the D = i case, in Sec. |IV[ the roots on the complex plane do not behave in the same way that they do in 
the four-dimensional case (Fig. 18); hence we cannot obtain the distribution of the roots by using the random number 
generating algorithm for the compactified M* pattern. 


VI. THE N DEPENDENCE AND THE LARGE -A LIMIT 


In this section we are going show what is the effect of changing the number of fields, i.e. iV, for particular dimensions. 
We can even compute the limit when A —>■ oo (cf. [T^[T3]) , which is equivalent to the spherical model [57]. To obtain 
the RG equation for large-A, we are going to rescale Eq. Q by A, and considering the new variables p —>■ p/N and 
u —>■ u/N. The derivative of the potential remains invariant under this rescaling u' —> = u'. As a first step, we 

divide the RG equation 0byA: 


dt 


u 

N 


{D-2) 




Ad 

1 + u' 


Ad 1 Ad 1 
W 1 + u' ^ W 1 4-A + 2pu" ■ 


(29) 


Next we perform the rescaling 


dtu = {D — 2)pu' — Du -b 


Ap 

1 + u' 


Ad 1 Ap 1 
~N~ 1 + u' ^ W 1 + u' + 2pu" ■ 


By taking the limit A —>■ oo the following terms remain: 


dtU = {D — 2)pu' — Du + 


Ap 

1 + u'' 


(30) 


(31) 


From this equation one can derive the VBF curves with the usual steps; just the A dependence is absent from the 
formula Eq. (14). 

In this section we cannot show of course the full A dependence; these are only a few checks for particular values of 
A. One should derive the root structure for each A, independently. However, these few examples could give us some 
idea what is going on when A is being changed. 


A. A dependence for 0(A) theories in D < 2 


Let us consider three different two-dimensional theories with field components A = 7 and A = 15 and A —> oo. 
The positions of their roots are shown in Fig. We can see in the two finite A cases that there seems to be a 
problem in the low orders of the expansion: the position of the roots will not satisfy the requirement of the M* 
pattern; however, as the order grows the pattern is restored. We need to mention that when a particular root gets 
very close to —1 in the next order sometimes it becomes complex, but its real part still stays around —1; hence one 
cannot see this accumulation exactly, but rather we can say this root ’’melted” into —1. This might be related to the 
fact that we already explained for the D = 3 case: the VBF polynomials get extremely flat close to — 1 because of the 
(1 -b m^)” factor in Eq. (14); hence it could mean some problem for the root finding algorithm to provide the right 
value. One way or another, a complex root does not define a physically sensible theory; thus we can still look at this 
pattern which will provide Eq. (20) as the probability density of the root positions when n —> oo. 
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(a) N = 7 


(b) N = 15 



FIG. 19: The N dependence of the root positions in D = 2. We can see that the M* pattern restores at larger order 
of the expansion for various N. Sometimes one will find complex roots close to —1 which can be considered 
unphysical. In this way it does not actually matter if we consider the limit of the distribution of the roots or it gets 
complex providing an unphysical situation. When the n —> oo limit is taken the only two roots that can be found 
are —1 and 0. In the large-IV limit no complex valued root has been found till the order n = 34. 


However, no complex roots occur for the largecase where we know that only the symmetric phase must be present, 
as it was proven analytically in in the framework of FRG. The fact that a root gets complex or accumulates at 
one of the stable fixed points —1,0 is irrelevant from the point of view of the physics: in both of the cases only — 1 
and 0 survive when n —>■ oo. For D = 1 we can find similar results. These were, of course, only a few checks on the 
N dependence and it should be considered case by case for every TV; however, from our experience we could assume 
that in n —>■ 00 roots inside the interval (—1,0) either get complex or they are arbitrarily close to —1 or 0 for any field 
component N. 


B. N dependence in 2 < D < 4 


In dimensions 2 < £> < 4 for finite N we will get a very similar result to the one which we obtained in Sec. |IV| For 
these values of the dimension, one has the richest fixed point structure of all; hence it should be checked case by case. 
However, here we will only focus on the integer valued dimension D = 3. For the large-iV in = 3 the positions of 
the VBF roots show a significant difference comparing them to the finite N (compare Fig. 20 and Fig. [^. We can 
identify the Wilson-Fisher fixed point just as we did it in Sec. EYl but in this case we cannot see the additional ” fake” 
fixed points at any truncation level. Also for the > 0 roots we can see differences: they start to be organised in 
the pattern which leads us to the distribution found in Sec. 
complex roots are found unlike for finite N. 

Comparing our results to mis], where in the large-A^ the fully analytic solution was derived without relying on 


HI Eq. (20). In this case it is also remarkable that no 
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the polynomial expansion, we can see that in the case of the analytic solution a line of fixed points starting from the 
Gaussian is found, whereas in our case, using the VBF, seemingly, we do not have that. However, we can observe 
a unique behaviour of the VBF curves, namely, that the multiplicity of the Gaussian fixed point grows with the 
truncation order. More precisely, we have 


n—2 

A„ oc (-l)"+i(m 2 )'^(")(l + m^r E 


(32) 


where "/{n) is an integer function of the truncation order n: 


7(u) 



if n = 2, 3 
if n = 4 
if n > 5. 


(33) 


This is the only case (that is found), which has the property of growing multiplicity of the root = 0. That 
might be a fingerprint of the line of fixed points which can be found in the full analytical solution presented in the 
aforementioned papers. 



5 10 15 20 25 30 35 

n 


FIG. 20: The roots of the VBF in = 3, iV —>■ oo are shown. Gomparing it to Fig. the most striking difference is 
that there are no ’’fake” roots between the root which is indicating the Wilson-Fisher (m^ = —0.388) and the one 
which is for the Gaussian, i.e. = 0. It is equally interesting that in this case we do not find any threshold beyond 

which the roots become complex in the > 0 region. 


C. N dependence in H > 4 

One of the most important results for theories D > A was the triviality (Sec. |^. Our finding is that the modification 
of the number of the field components N does not give any qualitatively different result for any such theory, when N 
kept finite. However, we can discover differences between the results when N < oo and V —?> oo for theories in H > 4. 
In four dimensions we will see the ” compactified” M* pattern of the real part of the roots as in Fig. The only 
question is whether the convergence on the complex plane towards the origin still holds. We present the results in 
Fig. for iV = 4, V = 15, IV = 40 and N —>■ oo. We can observe that the convergence to the origin slows down 
as N grows, however, for the N < oo cases, it is possible to fit a curve on them near the origin. The same form of 
the power-law is found for each N, that is, f{x) = ax^, just like in the case of V = 1 in Sec. The parameters 
are the following: a = 0.193 ± 0.001, b = 0.587 ± 0.002 for iV = 4; a = 0.249 ± 0.002, b = 0.586 ± 0.002 for N = 15; 
and a = 0.302 ± 0.006, b — 0.584 ± 0.004 for N = 40. It is interesting that close to the origin the exponent of the 
power-law seems to take the value 0.586 ± 0.002, and thus one can speculate whether it is universal. 

Regarding the V —>■ oo case, we did not reach the region where we could do the fit; however, the absolute value 
of the imaginary part reached its maximum and started to decrease. We expect the same behaviour (maybe with a 
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(a) iV = 4 (b) AT = 15 




Re 


Re rnP 


(c) iV = 40 (d) AT ^ oo 

FIG. 21: The convergence to the origin on the complex plane is shown for (a) A^ = 4, (b) = 15 and (c) N = 40. 

On (d) the A^ —>■ oo is presented; in this case the number of the data points is not enough to fit a curve on the points 
near the origin (the expansion order in this case is n = 53, and for the finite N cases n = 35). However, the same 
behaviour is expected: the convergence to the origin slows down as N grows. For the finite N cases a power-law fit 
can be found around the origin (the fit was performed on 24 data points), for which the exponent seems to be 

universal. For details see the text. 


different exponent) as for the finite N cases; however, to see that we would need to go to a higher expansion order 
than the one used in this case (n = 53). 

For theories in greater than four dimensions the fixed point structure does not change much compared to the results 
of Sec. |Vj only quantitative differences can be observed. We will present a few plots on the root positions in Fig. [22] 
Arbitrarily we have chosen the cases {D = 5, N = 6}, {D = 8, N = 10} and {D = 11, N = 12}. The roots are 
positioned in such a way that they accumulate more and more at the two stable roots —1 and 0, which is essential 
for triviality. There can be minor deviation from the pattern M*; however, at some point the pattern is restored 
completely, and one can obtain the triviality in the same way that we did it for the case D = 5, N = 1. 


A more interesting case that can be found is in the large-A^ limit. From Eq. (31) one can derive again the VBF 
curves in D > 4. In the present section we are going to study only the D = 5 theory for the large-A^ limit. In Fig. 
one can see the root positions for this model. The first thing that we can notice is that there is a stable line at 
= 0.139 which would signal a new fixed point solution. There has been a recent work [3S] on the topic of whether 
such a fixed point exists in models 1 < D < 6 and the answer was that, if there is such a hxed point it will provide 
an unphysical fixed point potential. In an earlier work |36j related to holography a similar conclusion was drawn. 

In the papers m an IR hxed point was found for an 0{N) symmetric theory with N + 1 scalars in = 6 — e for 
sufficiently large N, including a cubic interaction, too. It was argued that this IR hxed point of the cubic 0{N) 
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(a) D = 5, N = 6 


(b) D = 8, TV = 10 






• 



. ‘ 


• 

• 

♦ 

A • 



1 i 1 


t J 5 § 

* 1 ? 1 

: » ! 

• 

♦ 


■ A 

5 10 

15 

20 

25 


n 


(c) H = 11, TV = 12 


FIG. 22: The VBF roots of the 0(TV) model in (a) H = 5 and TV = 6, (b) Z1 = 8 and TV = 10, and (c) D = \1 and 
TV = 12. On each plot the accumulation of the roots at m? = — 1 and = 0 can be observed. For the roots 
w? < — 1 there can be small deviation from the M* pattern (most likely due to the extreme flatness of the VBFs 

around —1), but at some point the pattern restores completely. 


theory with TV + 1 scalars is equivalent to a UV fixed point of the 0(TV) model in the large-TV, which would imply an 
unexpected asymptotically safe behaviour of such theories. Now, since we have found such a fixed point, we need to 
check whether it gives a stable fixed point potential. Substituting the value w? = 0.139 into the VBF polynomials will 
provide the exact value for the corresponding coupling at that fixed point, hence defining the fixed point potential. 
The finding is the following: the highest order in the expansion is n = 46, and what we get is A„(to^ = 0.139) < 0 for 
2 < n < 46. This signals that the fixed point potential defined by this root is unbounded from below in agreement 
with the findings in . Since the RG equation of the effective potential can be considered to be exact in the large -TV, 
this unstable fixed point should be found also beyond LPA. In Fig.j^one can also observe that there is a gap between 
the roots corresponding to this newly found unstable fixed point, which can be filled in by considering the real part of 
the roots. In the same figure, considering the lower panel, we can see how the real parts of such roots behave. We can 
fit a power-law decaying curve again on the real part of the first root. This curve goes through the line rri^ = 0.139 
and in the asymptotic limit goes to zero. The imaginary part behaves pretty much in the same way as in the 0 = 4 
case (see Fig. [21] ). Here, we can use the same argument as for = 4, TV = 1: one needs to consider the positions of 
the roots on the complex plane in order to catch the underlying physics. According to our findings, the real parts of 
these roots exhibit the same M* pattern; hence for n —oo we can expect them to accumulate at = 0, leaving us 
with three fixed points, which are defined by the roots: = — 1 (since for < — 1 the root structure is the same 

as it was for the finite TV case), = 0 and = 0.139. Out of these three roots only the Gaussian {m? = 0) seems 
to provide a stable fixed point potential; thus triviality was found again. 





24 



FIG. 23: In the left panel the root structure provided by the real roots is shown. A stable line of the roots can be 
observed at the value = 0.139. In the right panel the same root structure is shown, but now the real parts of the 
complex roots have been taken into account, too. The real parts of the complex roots show the M* pattern, and 
also a power-law curve can be fitted on them: f{x) = ax^, with a = 0.807 ± 0.005 and b = —0.492 ± 0.002. In the 
n —)■ 00 we expect to find triviality again. For details see the text. 


VII. THE FRACTAL DIMENSIONS 
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(c) D — 3.1, N — 1 {d) D — 4.6, A —>■ oo 

FIG. 24: The root structure is shown for some fractal-dimensional cases. The usual fixed points are found for 4 < H; 
thus the most remarkable is (d), where we obtain a fixed point candidate for D = 4.6 in the large-A. However, it 
turned out to define an unstable fixed point potential just like it is shown in |35j . 
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We present a few examples of the fractal-dimensional cases in the intervals of the dimension that we have in¬ 
vestigated. In Fig. the results are shown for {D = 1.3,7V = 1}, {D = 2.6,7V = 1}, {D = 3.1,7V = 1} and 
{D = 4.6,7V —>■ oo}. In the first three cases the usual fixed points are found; however it could be possible to find 
additional fixed points for 2 < H < 4, but that would require a more detailed study of this interval of the dimension. 
For the two D values, which are in this interval, we found the WF fixed point, too, and it can be observed that they 
are situated at different values of for D = 2.6 rri^p = —0.34 and for D = 3.1 it oscillates around m^p = —0.16. 
The most remarkable case is when we set TV — )■ oo and choose D = 4.6 G (4, 6). Just as we saw it for the large-TV in 
V = 5, here, we find again a fixed point candidate at the value = 0.108. However, again in agreement with [35j . 
this fixed point is found to be an unphysical one, i.e. the root = 0.108 defines an unstable fixed point potential. 
(More precisely, it is a metastable potential.) 


VIII. CONCLUSION 


We investigated the fixed point structure of the 0(7V) model for various dimensions and filed components. In this 
paper we considered only the integer-dimensional cases in detail, but results for fractional dimensions are presented in 
the last section, too. For our analysis we used the LPA with a Taylor-expanded potential around zero VEV of the field. 
In this case from the fixed point equation one is able to express all the fixed point couplings through a polynomial in 
the mass of the theory; see Eq. (14). Each nth polynomial is obtained from the fixed point equation of the (n -|- l)th 
coupling; that is, the equation for 9(A„+i(A„, m?) = 0 can be solved for A„ = A„(m^), hence the name VBF (vanishing 
beta function) curve. A general property of these polynomials is that they always have a root at w? = — 1 and w? = 0 
corresponding to the convexity (or IR) and the Gaussian fixed point, respectively. These polynomials can be obtained 
for arbitrary components of the field TV and dimension D. To find a true fixed point at a given level of truncation, we 
established the following rule: let mg be a root of the polynomial A„(m^); then a physically well-defined (i.e. bounded 
from below) fixed point potential is given by the set of couplings {A„_i(mQ), A„_2(mQ),..., Ai(mg) = mg}, provided 
A„_i(mQ) > 0 [if it happens to be a root for the (n — l)th polynomial, too, then the same rule holds for the (n — 2)th 
polynomial, and so on]. Using this rule to find fixed points at a truncation level n gave us a nice opportunity to make 
considerations about theories in different dimensions. 

First, we considered theories for D < 2, and in particular, we analysed the case of 7A = 2 in detail. The finding is 
that for all the theories we considered up to and including dimension 2 they behave qualitatively in the same way 
for TV > 2. In such theories for continuous symmetries the Mermin-Wagner theorem must hold as it was shown in 
[2H] . We provided a statistical approach to show the validity of the theorem by considering the statistics of the root 
patterns. It was found that for such systems {D < 2) all the roots are in the closed interval [—1,0]. Since at every 
truncation level a new root appears, it seems to be paradoxical to prove the Mermin-Wagner theorem: we would not 
expect any root inside the interval but practically there appears (countably) infinitely many. A true fixed point in 
this interval signals a SSB phase; hence the MW theorem seemed to be violated. However, it turned out that this is 
only the case for finite n truncation. As n —> oo we can recover the MW theorem by simulating the position of the 

For discrete symmetry TV = 1 the 0(TV) model is 


roots and derive their distribution in the infinite limit; see Sec. HI 


equivalent to the Ising model which has a SSB phase \n D — 2. This result is also shown as a part of this section. 

In Sec. |IV| we investigated the theory for D = 3 with TV = 2, and in 2 < TA < 4 the finding is that the theories that 
have been considered have similar properties. We found an additional fixed point beside the Gaussian and the IR 
fixed point, namely, the Wilson-Fisher fixed point. Although these fixed points have been found, we could not clearly 
analyse the root structure for the > 0 region. It is also possible that one could find additional fixed points in 
2 < ZA < 4 for fractional dimensions, but it would require more detailed study. The main result of this section was 
that we were able to show the appearance of the Wilson-Fisher fixed point. 

For iA > 4 we found two qualitatively different results (see Sec. |^. However, essentially they lead us to the same 
physics, namely, to the triviality of the model; that is, only the Gaussian fixed point exists for such theories. For 
theories TA > 4 we analysed the D — 5 case, where the fixed point structure shows similarity to the D = 2 case, but 
now the roots only appear outside of the interval (—1,0); however, for the root pattern the statistical simulation can 
be used again if we compactify the real line of m^. In this way triviality can be shown. The situation for TA = 4 is 
different: here we need to consider the complex roots that appear during the calculations, but it can be shown that 
both the real and the imaginary part converges to the origin on the complex plane; hence leading to the triviality of 
the model. 

In Sec. |VI| we showed some results considering different values of TV. The most interesting outcome is connected to 
the recent results in [35l 1 ^ . where it was argued that in the 0(TV) model for TV —>■ oo if a non-trivial fixed point 
exists in dimensions 4 < TA < 6, then it defines an unstable fixed point potential. From our analysis we found such a 
fixed point, and indeed, it can be shown that it defines an unbounded effective potential. In the last section we gave 
some results for fractional dimensions, obtaining similar behaviour to the integer-dimensional cases. 
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According to our findings, we can draw the conclusion that although we used the local potential approximation during 
our computations, that seems to be enough to obtain the right qualitative physical results (Mermin-Wagner theorem, 
the presence of the Wilson-Fisher fixed point and triviality), when one extrapolates to infinite order the results 
obtained in a given order of the Taylor expansion. On the other hand, one has to be careful with such expansions, 
since as we saw, at finite level of the truncation it can generate fake hxed points from the physical point of view. 
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Appendix A: The proof of the nested formula 


We would like to prove the nested formula for the coupling constants in Eq. (13). The easiest to consider first is 
the case for the large-A^, D = 2, because the RG equation has the simplest form for these parameters [see Eq. iH)], 
however, it can be generalised to arbitrary N and D, which we will present after this simpler example. Without the 
loss generality, we set the constant Ajj = 1. Our starting point is the RG equation, and we start to differentiate it 
with respect to p: 


-2m(p) + w X , , , 
u'{p) + 1 

u"{p) 


dtu = 
dfu' = —2u'{p) — 

dtu” = —2u"[p) 4 
dtu”' = —2u^^\p) — 




K(p) + ir 

2u"(p)2 

{u'{p) + lf {u'{p) + lf^ 

6u"{p)^ Qu^^\p)u"{p) u^'^\p) 


{u'{p) + 1)" {u'{p) + 1)^ {u'{p) + 1)^ 


(Al) 


Let us rewrite the last three equation in the following form: 

dtu' = Fi{u') + g{u')u", 
dfu" = F 2 {u' ,u'') + g{u')u^'^\ 

dtu" = F^{u ,u\u") + g{u)u^'^\ (A2) 

where we dehned 


Fi{u')= -2m'( p), 
F2(m',m")=-2m"(p) + 


2u"{pf 


u^'^\p) 


{u'{p) + 1) (m'(p) -f 1) 
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We can hnd the following relations between Fi’s: 
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Of course, dFi/dp = 0 since in this case there is no explicit p dependence, but for the sake of consistency, we will 
indicate this term, as well. We can make the following statement for n > 1: 

= F„(u', u", u"',..., + g{u')u^'^+^\ 

+1^""' + - + (IJSt + 

(A5) 


We are going to show this by indu ction. Let us suppose Eq. (A51 is true. We are going to show that it holds for n +1, 
too. Let us differentiate Eq. (A5) once with respect to p. It yields 

dF, 




op ou' 

dp ^ du' ^ du" du' ) 

= E„+i(u', u", u '",...,«("+!)) + g{u')u^^+^\ 


(A6) 


where 
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/ 5F„ 9g(u') 




du' 


u" 


(A7) 


In this way it was shown that the rhs always depends on the highest derivative of u{p) linearly, e.g. in Eq. (A51. 

Now, if we look for the scaling solution then the Ihs vanishes; hence from Eq. (A5), 


9{u') 


Evaluating this expression at p = 0 gives 


•^n+l — 


E„(Ai, A 2 , A 3 ,..., A„) 
5(Ai) 


From this expression the nesting is straightforward: 

(-^15-^2; -^3 j • j -^n) -Fh(m ) 


'^n+l — 


5(Ai) 


g(w?) 


(A8) 


(A9) 


(AlO) 


Here, we used the notation m? = Ai. Thus, the formula of the nested couplings in Eq. (13) is proved for N ^ 00 and 
D = 2. 

Note that in the finite N case there are terms in the initial RG flow like pu' and pu". Evaluating the equation at p = 0 
is crucial to be able to neglect these terms, which would prevent us from performing the nesting. Thus, expanding 
into Taylor series around zero is the only case when we can define VBFs. 

In what follows, we will give a proof of the nested formula in Eq. (13) for arbitrary field components N and dimensions 
D. Again, we set the constant = 1. Let us start again by differentiating the RG equation [in this case Eq. 0 ] 
with respect to p: 


TV — 1 1 

a,u =(D- 2 )p.'(^) - Duip) + 
a,„' = {D- 2 )^"(rt + (D - 2)„'(rt - Du'{p) - 
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Let us express the equations above in the following way: 

dtu = Fi{u', pu") + g{u)u' + h{u', 2pu'){Zu' + 2pu^^^), 
dtu” = F 2 {u\u' , pu^^'*) + g{u)u^^^ + h{u ,2pu'){5u^^^ + 2pu^'^^), 


(A12) 


where 
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We can establish the relation again just like in Eq. (|A4[): 
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We can make the following statement for n > 1: 
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Let us suppose that Eq. (A15) is true for the nth term. Now, we will show that it is true for the (n + l)th term: 
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aM"+i 
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In this way, by induction, we could show that the statement of Eq. (A15) is true. Here (in the case of hnite N\ we 
have to set p = 0 and only then is it possible to do the nesting like in Eq. (AlO). From this point it is straightforward 
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to show this. 
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